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The rock-paper-scissors game is a model example of the on-going cyclic turnover typical of many 
ecosystems, ranging from the terrestrial and aquatic to the microbial. Here we explore the evolution 
of a rock-paper-scissors system where three species compete for space. The species are allowed to 
mutate and change the speed by which they invade one another. In the case when all species have 
similar mutation rates, we observe a perpetual arms race where no single species prevails. When 
only two species mutate, their aggressions increase indefinitely until the ecosystem collapses and 
only the non-mutating species survives. Finally we show that when only one species mutates, group 
selection removes individual predators with the fastest growth rates, causing the growth rate of the 
species to stabilize. We explain this group selection quantitatively. 

PACS numbers: 87.23.Kg, 05.10.Gg, 87.23. Cc, 87.18.Hf 



Introduction — When multiple individuals depend on 
a shared and limited resource, it is in the long-term in- 
terest of everyone that the resource is rationed to avoid 
depletion. However, in the short term it is in the in- 
terest of each individual to consume resources fast so 
as to gain a competitive advantage over their more pru- 
dent neighbors. Acting rationally to promote their own 
self-interest, each individual therefore increases consump- 
tion until the resource is depleted [TJ [5] . This dilemma, 
known as the tragedy of the commons, illustrates the 
need for restrictions on the use of limited resources to 
ensure sustainable development. The potential conse- 
quences of not regulating common property are particu- 
larly evident when it comes to issues such as overfishing 
and global warming. Surprisingly, even among primi- 
tive life-forms, such as bacteria and plants, prudent con- 
sumption of shared resources has been observed in large 
groups of competing individuals 3-6 . These individuals 
have no means of enforcing common restrictions, so the 
emergence of restraint must have an evolutionary origin. 

We here study a community of three species with a 
cyclic interaction where species 1 overgrows species 2, 
which overgrows species 3, which, in turn, overgrows 
species 1 (see fig. [T^,). Such intransitive systems, sim- 
ilar to the game rock-paper-scissors, have been identi- 
fied in numerous ecosystems, ranging from terrestrial and 
aquatic to microbial ecosystems [THT3"]. 

An interesting property of cyclic interaction is that 
growing fast will not help a species to gain biomass - 
it will help its predator! By growing too fast, the species 
will weaken the population of its prey thereby improving 
conditions for its predator (see fig. [Tjb-TIJ;) [THUS]. Thus, 
for the species as a whole it is advantageous to grow 
slowly, whereas each individual of the species will get a 
competitive advantage from growing fast. It has been 
observed, both in simulations and in experiments, that 
such a species will suffer the tragedy of the commons if 
they interact globally, but that biodiversity may be main- 




Figure 1: (Color online) a) We consider a system of three 
species with cyclic interactions. Species 1 overgrows species 
2, which overgrows species 3, which overgrows species 1. In- 
dividuals grow at different rates. The mean growth rate for 
species 1, 2, and 3 are denoted vi, V2, and V3. b) When 
species 1 grows faster than 2 and 3, species 3 will be more 
abundant than species 1 and 2. Here V2 = v$ = «i/5. c) 
When all species grow at the same mean rate, they will be 
equally abundant. Members of same species will self-organize 
to form clusters on the lattice. 

tained if the interactions are local [TBl - f2"2"] . The evolution 
of restraint must, therefore, depend on the spatial struc- 
ture of the system. However, a quantitative mechanism 
for this has yet to be identified. 

Model — In our model, each site of an L x L square 
lattice is occupied by a member from one of the species, 
1, 2, or 3. Initially, all individuals grow at the same rate. 
At each time step the following actions take place: 

• A random node i and one of its four neighbors j 
are selected. If i can overgrow j it does so with a 
probability Mj. Hereby, j becomes a member of i's 
species with the same growth rate Vj = 

• When i overgrows j, it might mutate by a small 
probability p m utate- Hereby, j will either become 
faster growing than i, Vj — (1 + j)Vi, or slower 
growing, Vj = +7), where 7 is a constant. 

Unless otherwise stated, we have used L = 200, 7 = 
0.02, and p m utate — 5 ■ 10~ 5 . This choice of parameters 
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ensured, that the range of growth rates within a species 
would only vary 2 — 4% at any given time. 

Results — Assume that the species grow at rates 
(vi,V2,v 3 ) and that the probabilities of finding each at 
any one lattice site are (j>i,P2, f>3)- In the mean field ap- 
proximation we achieve a steady state when v\pip 2 — 
V2P2P3 = v 3 p 3 pi , leading to O [33] 
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This steady state is stable for the spatially structured 
system. In the mean field model, where all sites are neigh- 
bors, the abundances of the three species will oscillate 
with larger and larger amplitudes until biodiversity is 
lost. It is, therefore, only possible to study the effect of 
mutations, when the system is spatially structured. 

If all three species are allowed to mutate, their growth 
rates will accelerate exponentially (see fig. [2^,). The 
growth rates for the three species stay approximately 
equal, so from ^ the species will remain equally abun- 
dant, as observed in fig. [2]3. 

When both species 1 and 2 are allowed to mutate, both 
their growth rates appear to accelerate linearly. Thus, 
from @, species 1 and 3 will continue to grow while 
species 2 decreases hyperbolically in size until it goes ex- 
ctint (sec fig. [2j>[2|i). After this, species 1 will quickly be 
overgrown by species 3, so only the non-mutating species 
is left. Due to the hyperbolic decrease in the abundance 
of species 2, the time before biodiversity is lost grows 
drastically with the size of the system. 

When only species 1 is allowed to mutate, it will evolve 
to a certain mean growth rate vi = (2.4 ± 0.1)«2. At the 
same time, species 3 will grow in size to p 3 = (2.6±0.1)pi, 
while the relative sizes of species 1 and 2 will remain 
about equal (see fig. |2b-|2t). The reason why V1/V2 differs 
from ps/pi is that the mean field approximation |2]) is not 
perfect. For larger system sizes the same mean growth 
rates are observed, but with a smaller variance. 

Interestingly, biodiversity is maintained in our system 
when either one or three species are allowed to mutate, 
but not if two species are. 

From ([2]) we get that the difference in growth rate 
between the fastest and the slowest species must stay 
bounded in order for the system to be stable. If not, the 
predator of the slowest species will decrease in size until 
it goes extinct. We therefore investigate the relative ac- 
celeration i ^ of the fastest species in the system. This 
is done by measuring how often a faster-growing mutant 
survives and, therefore, contributes to the total growth 
rate of the species. Naturally, this acceleration is propor- 
tional to the probability of mutation p mu tate , when this 
is small. More surprisingly, the acceleration is propor- 
tional to the square of 7 (see fig. [3^,). This is due to the 
fact that both the survival chances of a mutant and the 
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Figure 2: (Color online) a-b) When all species mutate at 
the same rate, their growth rates will increase exponentially. 
They will, therefore, stay equally abundant, c-d) When only 
species 1 and 2 are allowed to mutate, their growth rates will 
steadily accelerate relative to species 3. Thus, species 2 will 
become scarcer on the lattice, until it dies out quickly followed 
by species 1. Simulations are carried out on a 1000x1000 
lattice to better show the hyperbolic decline of species 2. e- 
f ) When only species 1 is allowed to mutate, it accelerates to 
a growth rate fluctuating around 2.4 times faster than species 
2 and 3. Consequently, species 3 grows to become 2.6 times 
more abundant than species 1 and 2. 



corresponding increase in growth rate are proportional 
to 7. The implications are that, for an ecological sys- 
tem where mutations of a broad range of magnitudes are 
expected to occur, the evolution of species will be domi- 
nated by big leaps. Small evolutionary improvements in 
fitness will most likely come to nothing. 

In fig. |3];, the relative acceleration of the fastest species 
is plotted as a function of the relative growth rates v\ jv 2 
and v\ /t>3 . These have been found by monitoring how of- 
ten fast mutants will succeed in outcompeting the species 
of growth rate V\, compared to slow mutants. Each point 
is an average over the introduction of a minimum of 80000 
mutants. It is seen that when the relative growth rates 
are large, the fastest species will decelerate. This is due 
to the better survival chances of the slow mutants de- 
creasing the growth rate of the species and promoting 
biodiversity. Multiplying all growth rates by a constant 
factor corresponds to changing the time scales, so the 
functional form of the figure is independent of the abso- 
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Figure 3: (Color online) a) The acceleration in growth rate is 
proportional to the 7 squared. Thus, the evolution of species 
in cyclic competition communities will be dominated by big 
leaps in growth rate, b) When all species are allowed to mu- 
tate, all species having equal growth rate is a stable fixed 
point of the dynamics. If species 1 is growing faster than 
2 and 3, faster mutants of species 1 will die out more fre- 
quently, thereby decreasing the growth rate, c) The relative 
acceleration in growth rate — ^S 1 as a function of the relative 

v\ at 

velocities v\jvi and V\jv%. When all species grow at the same 
rate, their growth rates will all accelerate. If a species is grow- 
ing at a much higher rate than the two others, slower mutants 
will have better survival chances than faster mutants, thereby 
decreasing the growth rate. Arrows indicate the flow of the 
system, when only one species is allowed to mutate. The thick 
line marks the theoretically optimal speed of species 1, given 
by equation |E]|. It is seen to agree well with observations for 
V2 « va- 



llate magnitudes of the growth rates. 

When all three species are allowed to mutate, their 
growth rates will remain equal to each other. If a species 
becomes faster than the others, its acceleration will de- 
crease, allowing the others to catch up (see fig. [3}d and 
[3j;). This constant increase in relative growth rate ex- 
plains the exponential acceleration in fig. [2^,. 

The case of two species mutating at the same rate cor- 
responds in fig. p3p to one species moving along the hor- 
izontal line Vi/v 2 — 1 and the other moving along the 
vertical line Vi/va — 1. Along both lines, the accelera- 
tion is positive, so the system will mutate to extinction 



as seen in fig. pb. If one of the mutating species becomes 
faster than the other, it will decrease its acceleration, 
thus allowing the other to catch up. 

Letting only one species mutate corresponds to moving 
along the line V\/v% = V\/vz in fig. [3b until the accelera- 
tion becomes zero at a value just below 2.5, as indicated 
by the arrows. This explains fig. [2^. 

Mechanism for deceleration — As demonstrated, 
the deceleration resulting from large relative growth rates 
is crucial for maintaining biodiversity. Using a simple, 
one-dimensional argument, we now derive a relation for 
the relative growth rates at steady state. 

Locally, a faster mutant will always have a competitive 
advantage over a slower mutant. However, mutants of 
different growth rates will have a tendency to separate 
spatially on the lattice [TT] . The faster mutants are then 
at risk of exhausting their neighborhood of prey, leaving 
it in an isolated cluster surrounded by predators (see fig. 
I2k-|4|;). To avoid this, the fastest species should allow 
time for its prey to grow through its predator, connecting 
it to a new cluster of prey (see fig. [4jzL-[4j") . At the optimal 
growth rate, a species will grow through a typical cluster 
of prey just when this becomes connected to a new cluster 
of prey. 

Since the abundances of the three species are given by 
Q, one can expect the typical cluster sizes to scale like 



(Ai, A 2 , A3) CX (y/vi, y^3, V^l)- 



(3) 



Imagine typical clusters of species 1, 2, 3, and 2 ar- 
ranged on a one dimensional line (see fig. |4g). As time 
passes, the two clusters of species 2 will overgrow species 
3 at rate v 2 . If species 1 is to overgrow the neighboring 
cluster of species 2 at the exact time when this becomes 
connected to the second cluster, we must have 
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When v 2 — U3, this has the solution v\jv 2 — v\jvz 
1 .'■'). which is in excellent agreement with the value of 
2.4±0.1 found in fig. [2^. In fig. [3fc the contour line of zero 
acceleration is seen to agree well with ^ when v 2 ps 1)3. 
The discrepancies for v 2 ^> W3 and v 2 <C v 3 arise because 
the cluster structures of species 2 and 3 disappear in these 
limits, which negates ^ allowing species 1 to survive. 

Discussion — Our results explain quantitatively how 
communities of primitive organisms, such as bacteria 
and plants, in cyclic competition can evolve to a state 
with moderate consumption of a limited resource. Even 
though individuals would get a competitive advantage by 
growing fast, groups of fast growing individuals locally 
deplete their prey. Since individuals growing at different 
rates have a tendency to separate spatially, this group 
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Figure 4: (Color online) a-c) Typical situation when species 
1 is five times faster than species 2 and 3. Species 1 invades 
a cluster of species 2 and quickly overgrows it all, leaving 
it surrounded by its predator, species 3. d-f) Typical situ- 
ation when species 1 is 2.4 times faster than species 2 and 
3. Before species 1 has overgrown all of species 2, this has 
grown through species 3 to connect to a new cluster of 2. 
g) Clusters of species 1, 2, 3, and 2 arranged on a line. In 
the space-time diagram species 1 grows with slope 1/vi and 
species 2 grows with slope l/v2- A typical cluster length of 
species 2 is A2 oc *Jvz and for species 3 it is A3 oc y/v\. The 
optimal speed for species 1 follows by demanding that a typ- 
ical cluster of species 2 is overgrown just when it connects to 
a new cluster of species 2. 



selection will favour moderate growth rates. Thus, the 
growth rate will be limited by the condition that clusters 
of prey should connect to other clusters of prey before 
being completely overgrown. The equation ([5| describ- 
ing this optimal growth rate is, to our knowledge, the 
first quantitative result in the field. 

The spatial structure of the system is crucial for the 
group selection. In a well-mixed system the species will 
start a race to extinction. In a locally structured system, 
such as a Petri dish or the ocean bed, biodiversity can 
be maintained if one or all of the species are allowed 
to mutate. If two species are allowed to mutate, they 
will increase their growth rate until the system becomes 
unstable. This interesting result has not previously been 
reported. 

If all species are allowed to mutate, but at different 
rates, the species with highest mutation rate will typi- 
cally become eliminated after a long transient dynamics 
where all species increase their growth rates enormously. 
In practice this drastic increase in growth rates will be 
limited by metabolic constraints. Therefore, also sys- 
tems with multiple evolving species should be stabilized 
before collapse. In conclusion, our results emphasize the 
importance of modesty in growth, as well as modesty 



in the ability to evolve towards larger growth rates. In 
particular, in the case of cyclic competition our results 
quantitatively explain how groups of primitive organisms 
may self-organize to a state of sustainable development, 
preventing the tragedy of the commons. 
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